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ABSTRACT 


Gravity data was collected from 82 seafloor and 41 land stations in 
a 334 sq km area between Santa Cruz and Ano Nuevo Point, California. 
The methods of data collection and reduction are discussed, with 
the introduction of a unique sequence pertaining to application of 
gravity data corrections. A complete Bouguer anomaly map is depicted 
and subsequently tied in with a previous survey of northern Monterey 
Bay. Isoline gradient analysis supports the concept that complete Bouguer 
anomaly profiles can be utilized to map granitic basement displacements. 
Complete Bouguer anomaly cross-sections are compared with correspond- 
ing profiles of seismic, well core, sea surface gravity, and magnetic 
data. Excellent correlation is exhibited between these profiles and the 
Palo Colorado-San Gregorio fault zone. Faulting within the Monterey Bay 
fault zone can be traced from analysis of related profiles. Supporting 
evidence of the purported dip-slip and strike-slip motion along the Palo 


Colorado-San Gregorio fault zone is presented. 
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L. INTRODUCTION 


Today, earth scientists can obtain crustal information directly from 
well drillings (8-10 km) or indirectly through seismic reflection and 
refraction profiles, magnetic anomaly, heatflow, or electrical measure- 
ments, and gravity measurements. The latter method is the subject of 
this paper. The measured factor, the acceleration of gravity, results 
from the gravitational attraction between the mass of the earth and a 
test mass. Changes in crustal density result in variations in the 
earth's gravitational field. The unit of acceleration, 1 emisece is 
called a gal (for Galileo) and it is possible to determine gravitational 
acceleration to within 0.01 mgal. 

This research is based on the occupation of 82 underwater ocean 
stations and 41 land stations during the spring of 1973. The bottom 
gravity data was collected through utilization of the Naval Postgraduate | 


School's research vessel R/V ACANIA (Fig. 1). 


A. OBJECTIVES 

The main objective’ of this study was to conduct a bottom gravity 
Survey in a fault zone area and produce gravity anomaly charts that 
could be tied in with a previous shallow water survey (Cronyn, 1973) and 
various land surveys (Bouguer gravity maps of California by Bishop and 
Chapman, 1967). Sea surface gravimetry in this region would prove difficult, 
if not impossible. The stability of the seafloor and minimum distance 
between the test mass (gravimeter) and the density contrasts of the 


crustal material lent high credibility to the results obtained. 
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No previous gravity anomaly charts of the survey area were available. 
The geological substructure of this part of Monterey Bay and adjoining 
continental shelf has hitherto been inferred from seismic profiling 
(Green, 1969 and 1970) and bottom samples (Martin, 1964). The present 
work was undertaken to verify, modify, and amplify the earlier geological 
analyses. It was hoped that the plotted gravity anomalies could aid in 
accurate delineation of the fault zones cutting through the area investi- 
gated. Thus, this research may assist scientists in making predictions 
of seismic activity. 

This report will first discuss the physical and geological setting 
of the surveyed area, prior to the development of the problem of gravi- 
metry itself. The experimental procedures are then discussed, followed 
by explanations of the numerous corrections necessarily applied in the 
reduction of the observed data. Finally, the gravity anomalies are 


presented in chart format, with accompanying analyses and interpretations. 


B. SURVEY AREA DESCRIPTION 

Santa Cruz on the northern edge of Monterey Bay, is located approxi- 
mately 104 km south~southeast of San Francisco. The survey area extends 
up the coast from senna Cruz to Ano Nuevo Point and from the 50-fathom 
curve to approximately one mile inland (Fig. 2). The area covered is 
about 334 Hs km. A gently sloping continental shelf is prevalent with 
small elevation gradients exhibited immediately along the seashore. An 
exception to this is the coastline between the town of Davenport and 


Ano Nuevo Point were cliffs are evident at the water's edge. 
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Figure 2. Survey Area Location 
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C. LOCAL GEOLOGICAL DESCRIPTION 

One of the most seismically active regions in North America, the San 
Andreas fault zone, lies just 24 km northeast of the surveyed area 
(Fig. 3). It probably represents part of a transform fault along which 
the Pacific and North American plates are slipping relative to each other. 
Scientists consider the entire coastal area west of the San Andreas to be 
geologically active. The region between the San Andreas fault zone and 
the Sur-Nacimiento fault zone to the southwest encompasses the Salinian 
block (Reed, 1933), which consists of continental crust dominated by 
granitic rocks, flanked on either side by oceanic crust of the Franciscan 
assemblage. North of Monterey Bay, faults in the Salinian block usually 
trend northwest-southeast and tend to offset the Cretaceous granitic 
basement rocks and overlying Tertiary strata (Jennings and Burnett, 1961). 
Offshore, the granitic basement of the Salinian block imparts a rigid 
block-faulting structural style to the overlying sediments (Hoskins 
andecrititiths, 19/1). 

The narrow Palo Colorado-San Gregorio fault zone (about 3 km wide) 
appears to come ashore at Ano Nuevo Point and ties in with the San 
Gregorio fault. Here the Miocene Monterey Formation appears to have 
been thrust to the southwest over the Pleistocene marine terrace deposits 
with a dip-slip component of about 3 m (Fig. 4). Four miles to the 
south, at Greyhound Rock, the Santa Cruz Mudstone and terrace deposits 
are vertically offset by 6 m at three separate but closely spaced faults 
(Griggs, 1973). Also, the terrace deposits on the sea cliffs where the 
San Gregorio fault goes out to sea are offset 2 to 3m (Clark, 1970). 
Offshore, the Palo Colorado-San Gregorio fault zone juxtaposes the 


Pliocene Purisma Formation and the upper Miocene Santa Cruz Mudstone. 
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Figure 4. Recent Fault Exposed in Sea Cliff at 
Afio Nuevo Point. (After Clark, 1970) 
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Fault plane solutions for five earthquakes since 1969 show nearly vertical 
fault planes and right-lateral strike-slip motion 6 m or less from the 
seatloon (Green et al., 1973). 

The Monterey Bay fault zone (10-15 km wide) also exhibits fault planes 
that are nearly vertical with right-lateral strike-slip displacement 
occurring on the northwest-southeast trending faults. Again the fault 
offsets often extend to within 6 m of the seafloor, generally displacing 
late Pliocene strata, some displacing Pleistocene deposits, and in a few 
cases, Holocene deposits. Thus they can be termed “geologically young." 
Most of these faults are downthrown on the landward side and are accom 
panied by drag folding (Green et al., 1973). 

On the continental shelf itself, the uppermost sediment is predom- 
inately sand with green mud in increasing evidence as Soquel Canyon is 
approached at the southeastern end of the survey area (Bishop and Chapman, 
1967). 

Visible proof of crustal motion in this region, historial accounts 
of earthquakes dating back to 1836, and the lack of reliable information 
as to the offshore geological substructure make this area very interes- 


ting to oOceanographers, seismologists, and geologists. 


D. PREVIOUS AREA INVESTIGATION 

Past geological and geophysical investigations of the area have been 
primarily reconnaissance studies (Shepard and Emery, 1941; Shepard, 1948; 
Martin, 1964; Curray, 1965, 1966; Rusnak, 1966; Martin and Emery, 1967; 
Hoskins and Griffiths, 1971; and Silver et al., 1971). The research was 
concerned with the general regional structure and major faulting as re- 


lated to the more extensive investigations of the central California 


18 





shelf. Martin (1964) utilized bathymetry and bottom samples as a basis 
for a geologic map of Monterey Bay but did not describe the area north 

of Santa Cruz. The land geology of the region has been thoroughly mapped 
by no less than 14 different persons starting with Johnson (1855), and 
most recently by Brabb (1970). The onland faults have been studied by 
Fairborn (1963), Sieck (1964), Durham (1965), Burford (1971), and Gilbert 
(1971). 

Clark (1970), and Evans and Lajoie (1971) have studied the northern 
end of the San Gregorio fault at Afio Nuevo Point. Hoskins and Griffiths 
(1971) suggested that faults across the Santa Cruz-Ano Nuevo Point con- 
tinental shelf did not affect rocks above a buried erosional unconformity 
of late Miocene age. Weber and Tinsley (1971) have also added insight 
as to the geological structure of the faults traversing Ano Nuevo Point. 

Jennings and Burnett (1961) and Jennings and Strand (1958) respec 
tively, compiled the San Francisco and Santa Cruz geological maps for 
the California Division of Mines and Geology. Bouguer gravity anomalies 
are plotted on the Santa Cruz sheet but those for the San Francisco 
Sheet have not at the time of this writing, been published. 

Bolt, Lomnitz, and’McEvilly (1968) postulated that the earthquakes 
occurring in the vicinity of Monterey and Santa Cruz in the 1800's and 
early 1900's may well have occurred in the two fault zones that cross the 
Survey area covered in this report. Griggs (1973) has chronolosically 
compiled the earthquake history of this region from 1836 to 1971. Most 
recently Greene et al., (1973) have depicted the crustal structure of 
the area in detail, paying close attention to the offshore faults in 


the Santa Cruz~Ano Nuevo Point sector. 


Ey 





E. BOTTOM GRAVIMETRY 

The measurement of gravity in the water covered regions of the earth 
is a difficult task. Gravity surveys at sea have been carried out from 
diving bells (Frowe, 1947), submarines (Meinesz, 1958), surface ships 
(Graf, 1958), and using a gravimeter enclosed within a waterprocf pressure 
chamber able to withstand many atmospheres of pressure at depth. Meter 
functions are remotely controlled from the attending surface ship through 
a number of separate electrical conductors located in an insulated, 
armored, oceanographic cable. A specialized winch with numerous slip- 
rings must be employed. 

Bottom topography plays an important role in bottom gravimetry in 
that the meter leveling system can operate only within a slope range of 
15°. The most important difference between bottom and sea surface surveys 
is the necessity of maintaining the research vessel above the meter and 
lowering enough cable so that the gravimeter will not be tipped over or 
dragged by the cable while the mass is in the unclamped mode. This could 
result in damage to the delicate precision spring suspension system and 
render the meter inoperable. Concurrently, navigation and accuracy of 
station plotting is ‘an additional inherent problem of sea gravimetry. 
Table I is a summary of the major differences between land and oceanic 


gravity surveys. 
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Table I. Major Differences Between Land and Oceanic Gravimetry. 
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If. SURVEY EQUIPMENT 


The first instrument used to measure gravity was a torsion balance 
(Cavendish, 1791); the pendulum was widely employed from about 1900 to 
the middle 1930's, and since then the gravimeter has been used almost exclu- 
sively. Because geophysical surveys are usually concerned with differences 
in gravity from one area to another, and not absolute values, there is 
limited need for an instrument that measures gravity directly. The 
gravimeter measures small variations in gravity, usually to a precision 


of about 0.01 mgal. 


A. LAND GRAVIMETER 

A LaCoste and Romberg Model G-08 geodetic land gravimeter (Fig. 5) 
provided by the U. S. Geological Survey (USGS) at Menlo Park, California, 
was used to conduct the coastal survey. Its characteristics and features 
are listed in Table 2. Figure 6 shows a simplified diagram of the meter's 
essential operating mechanisms. The zero-length spring exhibits a linear 
relationship between its elongations and compressions and increases and 
decreases in gravity itself (LaCoste, 1934). Rotation of the measuring 
screw brings the light beam to the equilibrium position by physically 
changing the location of the upper connection of the zero-length spring. 
Minimum meter drift is maintained by a thermister-transistor heater 


control system. 


B. UNDERWATER GRAVIMETER 
A LaCoste and Romberg Model H6G underwater gravimeter was used for 


the continental shelf portion of the survey. It was provided by the 


Jagd 








Figure 5. LaCoste and Romberg Model G-08 Geodetic Land Gravimeter. 
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MODEL G-08 GEO- MODEL H6G UNDER- 
DETIC LAND WATER GRAVIMETER 
GRAVIMETER 


Ger/2 in X 5a374 39 in triangular 
in X§ $e) 2 in base, 29 in height 


SIZE 
WEIGHT 
[a 


20 lb, including 350 1b plus addit- 
battery and carry- ional weight to 
ing case increase sink rate 


worldwide (7000 worldwide (7000 
mgal) mgal) 


OPTIMUM/REALISTIC 


+ 0.01/+ 0.04 mgal + 0.02/+ 0.10 mgal 
ACCURACY 





POWER 
SOURCE 


battery 15 vac, 60 Hz 


LEVELING manual micrometer automatic servo 


screws motors 


SYSTEM 





DRIFT less than 0.5 less than 1.0 
ngal/month mgal/month 





Table II. Land And Underwater Gravimeter Characteristics. 
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Figure 6. Simplified Diagram of the LaCoste 
and Romberg Gravimeter (after LaCoste, 
WMS GZ) 








Naval Oceanographic Office, Washington, D. C.. The underwater gravimeter 
is essentially a land meter mounted within a shell of two thick aluminum 
hemispheres. Its features, compared with those of the land gravimeter, 
are listed in Table 2. Again, internal temperature of the meter is 
maintained by a thermister-transistor circuit within the outer sphere. 
Clamping and unclamping of the gravimeter beam or weight, the meter level- 
ing system, and the measuring screw operation are all remotely controlled 
through a multi-conductor cable terminating at the control box on the 
research vessel. Figure 7 shows the complete underwater unit with the 
electrical termination located directly atop the pressure shell and the 
oceanographic armored cable protruding through the top of the termination 
housing. In Figure 8 the top hemisphere has been removed and the meter 
itself is visible, along with the leveling mechanisms and various other 
electrical components. 
1. Neste, Equipment 

Figure 9 shows a simplified diagram of the components of the 
auxiliary equipment. Primary power is provided by a gasoline engine 
which drives a hydraulic pump which in turn drives the cable winch and 
the A-frame through a -double set of two-way control valves. The hinged 
A-frame, Which provides a clear passage for the gravimeter over the ship's 
side, is shown in its fully-extended position with the meter ready for 
lowering in Figure 10. 

The permanent electrical termination of the cable is routed in- 
Side the hollow winch shaft and out through a slip-ring assembly which 
is connected via a cannon plug to the control box. The basic electrical 
power source is that of the shipboard 115 vac electrical system which is 
connected through an isolation transformer and a kepco rectifier to the 


control box. 
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2. Shipboard Installation 
The 33 m oceanographic research vessel R/V ACANIA was the plat- 


form from which the underwater survey was conducted. Its shallow draft 
of 3 m made it possible to take measurements in waters where larger ships, 
employing sea surface gravimeters, would dare not venture. 

The 4 tons of equipment was distributed about the after end of 
the upper deck, with the A-frame, its supporting plate, and the gravity 
meter located on the starboard side so that the person at the control 
box inside the dry lab could observe the raising and lowering of the 
meter (Fig. 11). The horizontal A-frame supporting plate was held in 
place with four 1/2-inch (1.27 m) bolts through the upper-level deck. 

The other equipment was affixed to the deck with lag bolts. Additionally, 
the reel frame was welded to the outboard fishplate via four solid steel 
bars. A modified hand brake for the real drum was installed as an added 


safety feature. 
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III. SURVEY PROCEDURES 


A. CALIBRATION OF GRAVIMETERS 

The land gravimeter used in the coastal survey features hardened 
micrometer screws and lever systems upon which the calibration factors 
depend. For this reason the calibration factors do not change perceptibly 
with time, eliminating the need for frequent checking. Nevertheless, after 
obtaining the meter at the USGS, a calibration run was made from the USGS 
headquarters in Menlo Park (USGS 1 JD) to Skeggs Point (USGS B-388). This 
route encompasses 5 benchmarks and a range of 137.20 mgal (Chapman, 1966). 
Subsequent reduction of the observed data on the IBM 360 computer at the 
Naval Postgraduate School yielded a difference of 137.13 mgal over the 
calibration range. 

Initial calibration of the underwater meter was carried out by LaCoste 
and Romberg technicians in Houston prior to shipment to Monterey. The 
excellent stable calibration characteristics previously mentioned for the 
land gravimeter are also inherent in the underwater model. Preceeding 
the author's survey, an additional calibration check was effected between 
the Woolard Airport Base WA-84 Station at the Monterey County Airport and 
Be chitark WH-29 at the end of the U. S. Coast Guard pier, Monterey Harbor 
(stations designated by Woollard and Rose, 1963). Slight modification 
to the readings due to recent construction at the airport station led to 
the estimation that the observed gravity difference between the airport 
and the pier stations was within 0.1 mgal of the value of 22.5 mgal as 


recorded by Brooks (1973). 
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B. COASTAL SURVEY 

A total of 41 shoreline stations were occupied during the survey. 

It was felt that a high station density along the shoreline would aid in 
the analysis of the data obtained from the offshore survey. A fairly 
constant horizontal range variation between stations was sought. Ultimate 
station selection was based on the accessibility of the USGS topographic 
map positions where elevation was recorded. Thirty stations were fixed 
at street or farm road intersections where the elevation was recorded, 
two were located at the maximum elevation of sea cliffs, and two were 
located at bridges where State Highway 1 crosses over creeks. An attempt 
was made to occupy as many USGS monumented benchmarks as possible. Of 
the seven benchmarks plotted on the maps, only three were found. The 
remaining four stations were located at the best estimated position of 
the plotted benchmarks, using surrounding topographical features and 
man-made structures indicated on the maps as references. The unoccupied 
benchmarks had most likely been destroyed by construction or were con- 
cealed by thick vegetation. 

It is estimated that the position accuracy of the coastal stations 
was at least twice as good as that of the underwater stations. Figure 
12 shows the plotted locations of the 41 land stations (letters) and 
Table 3 gives pertinent station location information. 

The first and last gravity readings each day were taken at the pier 


benchmark (WH-29). Observed meter drift was negligible throughout the 


coastal survey. 
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Figure 12, Land and Sea Station Locations 
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C. CONTINENTAL SHELF SURVEY 
The underwater survey encompassed a total of 82 sea stations. 
1. Station Selection 

Station locations were selected so as to cover the maximum amount 
of continental shelf in the survey area in the time available. A rela- 
tively constant spacing of approximately 2 km was maintained between 
stations. A grid of 100 stations was initially planned, but due to 
repeated equipment problems, rough seas, and limited ship availability, 
only the 82 stations were eventually occupied. However, a good represen- 
tative areal density was maintained and coupled with the extensive 
coastal survey, resulted in a reasonably large area being covered. The 
stations were numbered consecutively according to time of occupation 
(Fig. 12). The total 123 stations within the 334 sq km area yielded 
a station density of approximately one station per sq km, which was 
better than initially intended. Table 4 lists the underwater station 
location information. 

Stations at the eastern edge of the survey area were located as 
close as possible to the western edge of the shallow water area investi- 
gated by Cronyn (1972). It was felt this would simplify matters when 
attempting to tie in his pL anomaly isolines and this author's. In 
fact, two of Cronyn's stations were reoccupied for the purpose of com- 
parison. 

2. Navigation 

During this investigation, navigation was by visual bearings to 
prominent landmarks which were already accurately plotted on USGS maps. 
Generally, three bearings were taken at each station along with radar 


ranges to the coastline. Usually the estimated radius of accuracy of 
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each fix was small enough to be ultimately neglected during data reduction. 
A north-south plotting error of 0.2 km would result in only a 0.14 mgal 
difference in the value for theoretical gravity. Navisational and plot- 
ting procedures were carried out by two members of the ship's crew 
immediately upon lowering the gravimeter towards the seafloor. The ship's 
geographical position was kept approximately identical with that of the 
meter by heading the ship into the wind or prevailing swell prior to meter 
lowering and by maneuvering with the ship's engines thereafter. Sometimes 
cable had to be payed out to preclude the possibility of putting sufficient 
tension on it so as Co dra Or overturn the meter while it was on the 
bottom. 
3. Measurements 

Each day survey operations began with a measurement taken at the 
base station; ACANIA's two buoy harbor mooring next to WH-29, at the 
seaward end of the USCG pier. The absolute gravity measurement at WH-29 
was referenced to the counter readings obtained during several pier 
station occupations so that a formula could be employed to relate the 
daily base station readings with the known gravity value at WH-29. Upon 
reaching the desired location, the meter was lowered as fast as practicable 
to the seafloor. The meter lowering rate was usually a little over 
5 cm/sec. Bottom arrival was indicated when the depth counter units ceased 
fees nee) This value was recorded along with the fathometer reading. 
High speed leveling was initiated and flood and tilt checks made (meter 
leakage never occur and repositioning of the meter due to excessive bottom 
slope, was never necessary). After coarse leveling was completed the mass 
was unclamped. While fine leveling adjustments continued automatically, 


the beam position and gravity counter switches were manipulated 
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to null the beam position galyanometer needle. Hysteresis problems were 
avoided by always approaching the ultimate counter reading from low to 
high values. After obtaining and recording the correct counter value the 
mass was clamped by switching to the "deck" mode on the control box, the 
meter was raised and two-blocked, and the ACANIA was headed for the next 
assigned station. 

Under calm sea conditions and in intermediate depths (15 to 60 m) 
four stations could be occupied per hour. Shallow stations took longer 
because of galvanometer needle oscillations induced by the sea and swell. 
Deep stations proved to be time-consuming due to the large cable length 
required. 

4. Meteorological Ffforts 

It is felt there is a high reliability factor in the station 
plots since all survey operations took pjace during daylight hours. 
Although fog often prevails in the area during the spring months, no 
visibility problems were encountered throughout the survey. 

At wind speeds of 15 knots and greater more meter cable had to 
be let out. Even with constant maneuvering, the ACANIA could not be 
held in one spot. This problem resulted in a slight deterioration in 
position accuracy. The inevitable pitching and rolling often made 
operations more difficult and lengthy than predicted. On two separate 


occasions work had to be suspended due to high seas and/or strong swell. 
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IV. DATA REDUCTION 


In order to be useful to the geophysicist, observed gravity data 
must be corrected for station elevation, mass differences, the influence 
of nearby topography, and latitude. In general, all gravity values are 
reduced to a datum plane, which for this investigation will be taken as 
mean sea level (MSL). Due to the fact that absolute gravity values are 
rarely plotted on maps and charts, it is only necessary that anomalies 
be consistent as to the chosen datum. 

Much of the theory of the reduction of gravity data applies to both 
land and underwater surveys; but there are some important differences 
between the two. These will be explained in this section. Most of the 
actual numerical calculations were carried out on the Naval Postgraduate 


School's IBM~360 computer, programmed in Fortran langauge. 


A. OBSERVED GRAVITY 

As previously mentioned, a tie-in between WH-29 and ACANIA's harbor 
mooring was necessary to effect a working base station from which observed 
underwater gravity values could be derived through gravity difference 
calculations. The counter units recorded from the control box must be 
modified through a conversion factor in order to yield true milligal 
values. The conversion table provided by the manufacturers (LaCoste and 
Romberg, 1970) gives a single conversion factor of 1.03985 due to the 
fact that all counter values fell between 3300 and 3400 counter units. 

The absolute reference of 9/79891.7 mgal for WH-29 was correlated with 
the author's daily base readings (ba) taken at ACANIA's mooring. These 


values varied between 3323.67 and 3323.76 counter units. This variation 
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was due principally to earth and ocean tides and the difficulty involved 
in reading the meter scale during high sea and swell conditions. There- 
fore, the initial uncorrected observed gravity (C,) for the seafloor 


stations can be computed from the formula: 
G, = 979891.7 + (cv - ba) (1.03985) mgal, ‘e.) 


where cv is the control box counter value recorded at each station. The 
first half of each day's readings were referred to the morning pier reading 
and the afternoon readings were referenced to the evening counter value 
measured at WH=-29. 


For land stations the equation is: 
Gy = 979891.7 + (cv - bm) (cf) meal, (2) 


where bm is the pier benchmark reading taken with the land gravimeter. 
For the coast survey, the base station readings varied between 3405.47 
and 3405.67. Linear modification of the recorded counter values from 
the land meter was carried out through interpolation of conversion factors 
(cf in equation (2)) peculiar to meter G-08. All the calculations were 
done by computer using the Fortran program supplied by USGS. 

1. Earth Tide Correction 


The first correction that must be applied to G,, the earth tide 


1? 
correction (ET), is a result of the same forces that cause oceanic tides: 
the moon and the sun. Because the earth is not an infinitely rigid body, 
its inherent elasticity enables crustal deformations to occur due to the 
sun's and moon's gravitational attractive forces. These effects can vary 
the earth's radius as much as 390 cm in a few hours; yielding a net change 
of 0.1 mgal in gravity. Therefore, this correction can not usually be 
neglected in the determination of gravity anomalies. Again, a USGS com 


puter program was utilized to carry out the necessary calculations 
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for both land and sea stations, the input parameters being the geographical 
coordinates of the station. 
2. Instrument Drift Correction 

Readings were taken each day prior to sea survey operations and 
again upon returning to ACANIA's mooring. After removal of earth and 
ocean tidal variations, meter drift rate should be calculable between 
the times of base station occupation. Linearity in drift rate was not 
observed as the small amount of drift noted during the four months when 
operations were carried out was as much a positive factor as it was 
negative. Thus, meter drift was assumed to be negligible and was not 
considered in computing observed gravity. 

Drift for the land gravimeter was also assumed to be zero due 
to the characteristic calibration features inherent in the meter itself 
and the excellent correlation between previous USGS calibration runs and 
those of the author. 

3. Earth Curvature Correction 

The curvature correction (CC) is needed to compensate for the 
assumption of an infinite flat plate made in the ensuing Bouguer correction. 
Since terrain corrections were carried out to distances of 160 km from 
each station the flat earth assumption would ultimately result in signi- 
ficant errors. The USGS formula used was: | 
CC = ~ 1.376 x 10 *(z-2,) + 3.049 x 107” (2-2) = ee Bi aye 10°17 (2-2,)°, (3) 
where z is gravimeter depth and Ze tidal height, both in meters. Curva- 
ture corrections for the sea stations varied from -0.02 mgal at a 12 m 
depth to ~0.13 mgal at 90 m. Corrections for land stations were made 


from the equation: 


4 Oia =I 


CC = 1.376 x 10 H ~ 3.049 x 10 “H’ + 1.110 x 10 ~°H meal, (4) 
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where H is station elevation in meters above MSL. A maximum correction 
Gof 0.16 mgal corresponded to a 110 m elevation. 
It is now possible to calculate the observed corrected gravity 
(OG) from the equation: 
OG = G, ap OWE 8B e OOyr (3) 


where D is the meter drift (taken as zero in this work). 


B. THEORETICAL GRAVITY 

For the purpose of geophysical surveys, a reference ellipsoid has 
been adopted for use in calculation of the theoretical gravity (THG). 
This is the value of gravity expected if the earth were an ellipsoid of 
revolution fitted as closely as possible to MSL. 

Numerous variations in the equation for theoretical gravity as given 
by the International Gravity Formula have been derived. The equation and 
coefficents used by the author were those based on the formula of Heiskanen 
and Vening Meinesz (1958): 

THG = 978049.0(1 + 0.0052884 stall a0. 0000052 TS mgal, (6) 

Where L is the latitude. From differentiation of equation (6) one obtains 
an average south-to-north increase (in the Northern Hemisphere) in THG of 
0.81 mgal/km. This variation is due to both the decrease of the centripetal 
acceleration, due to the earth's diurnal rotation, from the equator to the 
poles and to the difference between the polar and equatorial radii. This 
variation in gravity with latitude constitutes the so-called "latitude 


correction." 


C. ADDITIONAL GRAVITY CORRECTIONS 

Four additional modifications have to be made to the underwater 
gravity data so as to obtain a complete Bouguer anomaly (CBA). In order 
of increasing difficulty they are: the free~air correction, the initial 
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and secondary Bouguer corrections, and the terrain correction. The 
most logical sequence for understanding the physical significance of each 
of the applied corrections is that prescribed by Andrews (1973), which 
is followed henceforth. 

1. Initial Bouguer Correction 

This step essentially removes the gravitational effect of the 

water above the meter, a problem peculiar to underwater gravimetry. An 
infinite "Bouguer plate" of water is assumed to lie over the meter, with 
properties of density c and thickness ‘zl. The equation relating these 


parameters to the correction (be, ) is: 


be, = 270 Gz (7) 


; ; ; =§ 
where G is the universal gravitational constant (6.670 x 10 
3 2 ; é ean 
em /g-sec ). This correction will be positive because the water above 
the occupied station exerts an upward attraction on the gravimeter. 


Given a water density of nas O02 7 oterie (this value is used throughout 


this paper) and measuring the depth z in meters, equation (7) becomes: 


be, = 0.043042 meal. (8) 


As previously stated,’ this correction does not apply to land stations 
as the. surrounding medium is air of negligible density. 
2. Free-Air Correction 

Next, the gravity station must be repositioned to the approxi- 
mated reference ellipsoid, that is MSL. If a gravity station were 
located exactly at MSL, this correction would not be required, but in 
this survey all stations were either below or above MSL. At MSL the 
free-air gravity gradient is 0.3083 mgal/m (Heiskanen, 1967). The 


general formula for the free-air correction (FAC) for underwater 
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stations is: 


FAC = - 2GM(z - z,)/R° A (9) 


WieLewielSmiEnercarth Simass (5.9/ x TO 1a). and R the earth's mean radius 


“6.571. x 11) a) (MacDonald. 1966). Equation (9) reduces then to the form: 


FAC = - 0.3083(z - Zz.) meal, (10) 
where z and Ze are measured in meters. 
Since station depth and tide level are the only variable 
parameters in equation (10) it is clear that the accuracy with which 
they are measured is critical. A transducer inside the bottom of the 
gravimeter sphere measures differences between surface and seafloor 
pressure. These differences are indicated by a counter in the control 
box which is related to depth by calibration tables provided by the meter 
manufacturer. 
The tidal heights were computed using the average between local 
stations at Santa Cruz and Ano Nuevo Island, related to San Francisco, 
the reference station for both. Since all tide tables use mean lower 
low water (MLLW) for their datum plane, it was necessary to determine the 
vertical distance between MLLW and MSL. This difference was found to be 
0.884 m (Coast and Geodetic Survey chart 5403). 
For underwater gravity stations, the FAC derived from equation 
(10) is negative since the meter is always below MSL. However, this 
does not hold true for land measurements, for which the following equation 
is used: 


FAC = 0.3083H meal, (11) 


where H is the station elevation in meters. In this case the meter is 


always above MSL and hence the FAC is positive. 
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3. Secondary Bouguer Correction 


This correction involves filling in the region between the 
reference spheroid (MSL) and the station elevation with rock of appro- 
priate density. In this way seafloor gravity observations can be corre- 
lated with land observations. The formula for this correction (be,) for 


Bo cCabloom Station is: 


be, = 2n 0 Gz - Zs (12) 


where om is the mean density of crustal rock in the Bouguer plate (Dobrin, 
1960). Numerical substitution for G, and using a value of 2.6/7 ale for 


C. (a value to be used throughout this paper) yields: 


be, = Q0.1119(z - Zz.) mgal, (13) 


where z and Zz are again measured in meters. This correction fills the 
space from the actual station depth to MSL with a uniform infinite plate 


of mean crustal density. For underwater stations, bc, is a positive 


2 
correction as mass is being added below the reference ellipsoid. 


For land stations the following equation is given: 


be, = - 0.1119H meal, (14) 


where again H is Beacon elevation in meters. This correction will be 
negative because we are effectively removing the material between MSL 
and the eomraelevat tone 

Combining the initial and secondary Bouguer corrections yields 
the complete Bouguer correction (BC). For the underwater stations the 


combination of equations (8) and (13) gives: 


BC = 0.04304z + 0.1119(z - z) mgal. (15) 
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For the land stations equation (14) is the complete Bouguer correction. 
Figure 13 schematically portrays the procedures for application of the 
free-air and Bouguer adjustments. 
4, Terrain erection 

By far the most tedious part of gravity data reduction is that 
of calculating the terrain correction (TC). The Bouguer correction 
assumes that the topography surrounding the station is that of an infinite 
flat plate. The terrain correction evaluates the error in the Bouguer 
correction due to undulations of the terrain about the plane through the 
station. For land stations the Bouguer correction overestimates the 
gravitational attraction of the mass below station elevation because 
it ignores voids in this space. For underwater stations, the attraction 
is underestimated since the Bouguer plate only extends down to the original 
Station depth, neglecting any voids below. The importance connected with 
this portion of the terrain correction in the author's survey area can be 
considered a direct function of the close proximity of large-scale 
bathymetric features. These include the Monterey Submarine Canyon and 
deep Pacific Ocean abyssal plains to the south and southwest, respectively. 
Therefore, the corrected gravity value obtained through application of 
the Bouguer correction will in actuality be too low (too much was subtracted 
for the land stations and too little was added for the seafloor stations). 
Furthermore, the upward component of the gravitational attraction of the 
mass above the plane through the station, which tends to lower observed 
pravity, 15 1lenored in the Bouguer correction. Again, it can be seen 
that this is an important factor in reduction of the author's data due 
to the proximity of the Diablo, Santa Lucia, and Santa Cruz Mountain 


Ranges. In summary then, topographic elevations above and depressions 
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below station location both act in the same sense, i.e., to reduce the 
gravity reading, and the gravity values modified only by the Bouguer 
correction will always be too low. Corrections for the effect of terrain 
elevation variations with respect to the plane through the gravity 
station will always be positive. 

The topographic attraction is most conveniently estimated by 
dividing the area around the station into segments. The terrain effect 
for any segment is then a function of the difference, whether positive 
or negative, between the station elevation (or depth) and the mean 
elevation (or depth) of the segment. Segmentation of the area surrounding 
each station was effected through the uSe of transparent templates 
(graticules) made up by the author. Three different maps were employed, 
requiring three separately scaled templates (1:24,000; 1:210,668; and 
1:820,000). Each template, when centered on a station, divides the 
surrounding terrain into compartments formed by radial lines through the 
Station, intersecting concentric circles centered at the template axis. 
The lettered areas between concentric circles are 'zones'' which are 
sectioned into numbered "compartments" by the radial lines passing through 
the station. Zone A, with two compartments, has an outer radius of only 
2m, while Zone 0, the most distant, has 28 compartments and an inner 
radius of 98.9 km and an outer radius of 166.7 km. After positioning the 
template over the desired station, the average elevation (or depth) of 
each compartment above (or below) MSL is visually estimated. The absolute 
numerical difference between station elevation and each of these compart-— 
mental values is then tabulated. With these values appropriate tables 
are entered to give the additional vertical gravitational attraction at 


the station for the compartments. Each compartmental correction is then 
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summed to give the zonal corrections which in turn are totalled to give 
the entire terrain correction for the station in question. There are 
199 compartmental corrections to tally for each station. 

The terrain correction tables used for this survey were based on 
those of Hayford and Bowie (1912). Bullard (1936) modified the Hayford- 
Bowie tables by enlarging the more distant compartments and reducing the 
compartments in the zones near the station. This eliminated some of the 
labor involved in use of the tables. In turn, the USGS modified Swick's 
(1942) work which was based on Bullard's tables. These then are the 
tables actually used in this research. A crustal density of 2.67 iene 
was assumed and a 0.615 multiplication constant was applied to oceanic 
compartments in which the average bottom depth was greater than that of 


the station. This proportionality factor was derived from the equation: 


eos {ey 


cw _ 2.67 — 1.027 _ 
ar 0.615 . (16) 


Computation of the corrections necessary for the land stations 
is relatively straightforward. However, care must be taken when averag- 
ing depths in the oceanic compartments. Use of the 0.615 factor for the 
entire compartmental column from land station elevation to seafloor fails 
to take into account the existence of air between the station and MSL. 
Since the eerie vertical distance from land station elevation to MSL 
was small in comparison to the overall compartmental depths, simple 
multiplication of the tabular values by 0.615 proved to be sufficiently 
accurate for all oceanic segments when calculating land station terrain 


4 


corrections. 
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Precise calculation of terrain corrections for underwater stations 
requires a two-step procedure. Because seawater surrounds the station 
instead of air, and because some of the nearby topography is above MSL 
due to the close proximity of the coastline, it again becomes necessary 
to apply a weighting factor to the terrain correction tables. 

Step one involves filling in the voids below station depth with 
rock of density 2.67 Pict (area 1, Fig. 14). Since water of density 
£027 on cm already fills these areas, the factor of 0.615 must be applied 
to the tabular values extracted for each ocean compartment with a bottom 
depth greater than that of the station itself. 

The next step entails removing the effect of the crustal rock 
lying above the meter (areas 2 and 3, Fig. 14). This is done in two 
operations since the rock in area 2 is immersed in seawater, while that 
in area 3 is surrounded by air. First, the total correction for areas 
2 and 3 is calculated assuming an average crustal density of 2.67 ak. 
As usual, this part is added to the observed gravity, as it represents 
mass removed from above the meter location. Figure 14 illustrates that 
at this point too much mass has in fact been removed since in area 2 a 
fraction of the terrain effect of the rock is counteracted by the overlying 
seawater. Therefore, part two involves computation of the effect for 
area 2 assuming a density for seawater of 1.027 Sak The total terrain 
correction for area 2 must be multiplied by the factor 0.385 derived from 
the equation: 


Dale 027 


wx 0.385 . (17) 


6} 

Sul 
Oo 

rc 


This second part of the correction should be subtracted from the observed 


gravity, as water of depth z (area 2) is added above the meter. 
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Figure 14. Schematic Diagram Showing Areas Involved in Terrain Corrections. 
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In summary, the total terrain correction for each seafloor station 
involves the following: 

(1) fill in area 1 with rock instead of water, 

(2) tit loam garea 2 with waterpinstead of rock, and 

(3) remove the effect of the crustal rock above MSL in area 3. 
In practice, the station depths (z) were small in comparison to the 
elevations of the terrestrial features (area 3); hence, the contribution 
of equation (17) was negligible. 

The possible error involved in calculating terrain corrections 
is estimated to be approximately + 0.02 mgal/zone which could ultimately 
result in a maximum possible error of + 0.35 mgal for the entire station 
correction. 

The maximum terrain correction was 6.53 mgal for Station LL located 
almost 2 km inland at the base of a steep-sided valley, while the minimum 
correction was 2.36 mgal for Station A located in downtown Santa Cruz. 

Sea-surface gravimetry usually yields massive banks of data com- 
pared to the relatively small amounts of bottom gravity data available. 
Correlation of the two is often a valuable tool to the geophysicist. 

With this in mind it is probably more realistic to first calculate the 
terrain correction for the bottom data prior to computation of the mass- 
adjusted free-air anomaly (to be discussed in the next section). This 
procedural alteration is proposed by Andrews (1973), based on the concept 
of downward continuation (Peters, 1949; Trejo, 1954) which asserts that 
the terrain correction for bottom gravimetry is greater than that for 
sea-surface gravimetry at a finite distance above the same location. 

Tables V and VI list observed and theoretical gravity and 
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D. GRAVITY ANOMALIES 


The difference between observed corrected gravity and theoretical 
gravity is an anomaly. When isolines are drawn to the gravity anomaly 
values plotted at their corresponding stations, analysis may reveal 
local and regional trends, which reflect geological substructure and 
density variations. Four different types of gravity anomalies were 
calculated for this survey. 

1. Free-Air Anomaly 

The vertical gradient of gravity can be approximated as the basic 
cause for gravity differences between adjacent stations. To a close 
approximation the vertical gradient is linear and is independent of 
latitude (see Equations 10 and 11). Application of the free-air correc- 
tion (FAC) (previously derived) results in the corresponding first-order 
anomaly known as the free-air anomaly (FAA). It is called "free-air" 
because the theoretical anomaly is computed as if the gravity measurement 
were made at MSL but without taking into account the attraction of the 
crustal material (or seawater) between the actual topographic elevation 
at the station and MSL; that is, as if the gravimeter were suspended free 


in the air. This anomaly is given by the equation: 
FAA = OG + FAC —- THG. (18) 


The FAC is always negative for bottom stations and positive for land 
stations above MSL. 
2. Simple Bouguer Anomaly 
When the Bouguer correction is needed ee the free-air anomaly, 
the simple Bouguer anomaly (SBA) results. The SBA is of little use in 


areas of prominently rugged terrain, but where uniform topography prevails 
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(e.g., Gulf of Mexico continental shelf) this anomaly can be utilized 
for gravity survey correlations and interpretations. The simple Bouguer 
anomaly is calculated as if the material under the station was of infinite 


horizontal extent. It is given by the expression: 


SBA = 0G + FAC + BC - THG, (19) 
Ox: 

SBA = FAA + BC. (20) 
For land stations BC (BC = bcos Equation 14) is negative; for seafloor 


stations BC (Equation 15) is positive. 
3. Complete Bouguer Anomaly 

Upon application of the terrain correction (TC) to the simple 
Bouguer anomaly the complete Bouguer anomaly (CBA) alee This anomaly 
is usually used to make composite maps and charts of isolines to tie in 
separately surveyed areas. The plotted CBA isolines should reveal lateral 
crustal density variations and near-surface structural non-conformities, 
as well as variations in the depth to the Mohorovicic discontinuity. Due 
to the relatively small area involved in this survey, changes in MOHO 
depth were considered to be insignificant. ae ee of gravity data 


reduction is given by: 
CBA = 0G + FAC + BC + TC - THG, (21) 
ors 


CBA = SBA + [Cane (22) 


4. Mass-Adjusted Free-Air Anomaly 


In order to be able to correlate bottom gravity measurements 


with sea-surface values, the mass-adjusted free-air anomaly (MFAA) is 
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introduced. Three steps are required: 
(1) remove the upward attraction of the overlying 
seawater (be,); 
(2) adjust the gravity as determined from the meter from 
its bottom value to that on the reference ellipsoid (FAC), and 
(3) fill the ocean (mow air below the meter repositioned at 
the reference ellipsoid) with seawater using the following 


mass-adjusted free-air correction: 


MFAC 


li 


ang Gz ~ z.)- 23) 
Numerical substitution for G and g yields: 


MFAC 


0.04304(z - Zz.) ngal, (24) 


Where z and Zz, are measured in meters. This correction will be positive 


as the underlying mass of water will increase the downward gravitational 
attraction when the meter is relocated at MSL. Thus, the MFAA is given 
by: 

MFAA = OG + be, + FAG Se MPAC = THG. (25) 


As previously mentioned for underwater stations, bc, is positive and 


a 
the FAC is negative. Obviously, this anomaly is peculiar to bottom 


gravimetry and has no terrestrial counterpart. 


In summary then, the forementioned corrections will have the 


signs: 
LAND STATIONS UNDERWATER STATIONS 
BC ~ es 
FAC + - 
TC 5 + 


Gravity anomalies for both land and sea stations surveyed are 


tabulated in Tables VII and VIII. 
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V. DATA PRESENTATION AND DISCUSSION 
A. GENERAL 


The computed values of the complete Bouguer anomalies for both land 
and sea stations were plotted at their respective positions on a compo- 
site large-scale USGS chart. IJIsolines were drawn by hand utilizing a 
2-mgal spacing. The final result shown in this report exhibits the most 
logical fit between the data and previous knowledge of the regional and 
local substructure. 

Inaccessability to precise basement depth data in the survey area 
precluded prediction of an accurate regional trend. Superficial exam- 
ination based on a speculative basement depth prediction indicates no 
Significant linear regional trend. Thus, no attempt was made to develop 
two-dimensional profile models of the area. 

The values calculated for the mass-adjusted free-air anomalies for 
all the underwater stations were also plotted on the same charts used 
for CBA analysis. Isolines were manually drafted, again at 2-mgal 
intervals. Both isoline charts were tied in with the previous gravity 


survey values for northern Monterey Bay calculated by Cronyn (1973). 


B. ERROR ANALYSIS 

Table IX lists what the author believes to be the maximum possible 
errors inherent in calculation of the complete Bouguer anomaly. It is 
felt that ultimate CBA values are accurate to approximately + 1.0 mgal 


for the sea stations, about five times more accurate than results from 
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sea surface gravimetry. As expected, it was possible to compute the 
CBA for the land stations with far greater accuracy, It is estimated 


Eiake they are@accurate to within + 0.55 meal. 


C. COMPLETE BOUGUER ANOMALY DISTRIBUTION 

Although plotted gravity anomalies reflect, among other things, 
variations in depth to the Mohorovicic discontinuity, it is believed 
that the small areal extent of the survey excludes the possibility of the 
existence of large slopes at the MOHO. Therefore, the assumption is made 
that the CBA reflects crustal deformations and general depth to the 
granitic basement. 

Figure 15 is a small-scale plot of the CBA isoline distribution for 
the survey area. The gravity high in the center of the area lies well 
within the Palo Colorado-San Gregorio fault zone. The high may be in- 
terpreted as an area of shallower depth to the basement complex Santa 
Lucia eepailtacy iin respect to the surrounding portion of the continental © 
shelf as will be shown later by seismic reflection profiles. 

The anomalous low centered off Davenport is situated on the eastern 
side of the fault zone and probably signifies a down-thrust portion of 
the zone, that is, the granitic basement is deeper there than in the 
ceo un dine areas. 

Northward, the large isoline gradients again reflect the presence of 
the Palo Colorado-San Gregorio fault zone extending into Ano Nuevo Point. 
The fact that the isolines do not run parallel to the direction of the 
suggested fault zone at first seems misleading. However, reports of 
eravity investigations inland near Ano Nuevo Point indicate similar CBA 


PSoline EeendsmGimmemm@larcs USGS oral commun.s,91973) . 
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Figure 15. CBA Distribution for the Continental Shelf and Adjacent 
Coastline Between Santa Cruz and Ano Nuevo Point. 
(Values in milligals) 
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The trending of the onshore isoline characteristics is based entirely 
on the author's data but the high station density of the coastal survey 
portion of this report lends much credibility to the CBA distribution on 
land (Fig. 15). 

The intensive curvature of the isolines to the north and south of the 
centralized gravity high may support considerations of dip-slip fault 
motion (and also possible strike-slip movement) along the Palo Colorado- 
San Gregorio fault zone. 

The tight gradient of the isolines in the southeastern part of the 
survey area is in excellent correlation with the Monterey Bay fault zone 
(Greene, 1973) where seismic profiling has indicated faults extending 
up to the base of the thin Holocene depositional layer. 

The large positive CBA values near Santa Cruz coincide with the out- 
cropping of the granitic basement complex north of Naturai Bridges State 
Beach. Also in evidence is the north-northeast, south-southwest slope 
of the basement offshore. The isolines, the trend of which is confirmed 
by the Santa Cruz outcrop, are perpendicular to the proposed downslope. 

A tie-in with the Bouguer Gravity Map of California was only possible 
in the land region near Santa Cruz as the Santa Cruz sheet does not extend 
beyond 37° N latitude. The general location of individual CBA isolines 
of the Santa Cruz sheet agreed well with the coastal and shelf stations 
of the author. However, once offshore the small extensions of the iso- 
lines on the Santa Cruz sheet do not conform exactly to the curvature 
of the author's isolines. Undoubtedly, this disagreement results from 


the previous land data being extrapolated to an extent beyond the shore- 


line. 
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D. MASS-ADJUSTED FREE-AIR ANOMALY DISTRIBUTION 

Figure 16 depicts the mass-adjusted free-air anomaly chart for the 
Santa Cruz-Afio Nuevo Point continental shelf. As expected, the general 
isoline trends run in the same directions as those of the CBA chart with 
the exception of the disappearance of the iow previously encountered in 
the CBA distribution. This is undoubtedly due to the additional numerical 
modifications involved in the calculation of the three-step mass~adjusted 
free-air anomaly introduced earlier as equation (25). 

The isolines terminate at the shoreline because there is no criteria 
for calculation of the MFAA at land stations. The MFAA plot is presented 
in this report solely for the purpose of possible correlation in the 


future, with sea surface gravity data. 


\ 


E. NORTHERN MONTEREY BAY TIE-IN 

Utilizing a pantograph to effect a 1/8th reduction in scale, composite 
charts were drawn using this author's data along with that of Cronyn's 
(1973) northern Monterey Bay survey. 

1. Complete Bouguer Anomaly 

Figure 1/7 portrays the composite CBA distribution resulting from 

connection of the northern Monterey Bay gravity values with those computed 
by the author. Isoline tie-ins were easily accomplished for the greater 
magnitude CBA values (10 through 25 mgal) in the northern part of Cronyn's 
area. At the extreme southeastern edge of the author's area correlation 
was not precise, although the anomalous trend was easily maintained and 
the extension of the author's isoline values to the east was readily 


carried out. 
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Figure 16. MFAA Distribution for the Continental Shelf Between Santa 
Cruz and Ano Nuevo Point. (Values in milligals) 
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The gravity low sited approximately 10 km south of Santa Cruz 
(Fig. 15) is substantiated in Cronyn's analysis and appears to be an 
area of downdropping of the basement complex. This coincides with faulting 
that extends to the base of the Holocene deposits within the Monterey Bay 
fault zone (Greene, 1973). 

The abrupt north-northwest, south-southeast trend of all but the 
10-mgal isoline to the south of Terrace Point is probably indicative of 
the previously postulated strike-slip features within the Monterey Bay 
fault zone. 

2. Mass~Adjusted Free-Air Anomaly 

Figure 18 represents the overall mass~-adjusted free-air anomaly 
distribution for northern Monterey Bay extending northwestward on the 
continental shelf through the author's survey area. Tie-in with Cronyn's 
data was easily accomplisned. This composite chart has been produced for 
the express purpose of comparison with and addition to, sea surface gravity 


data. 


F. CROSS-SECTION ANALYSES OF GEOLOGIC SUBSTRUCTURE 

Data relating to subsurface structure, other than that obtained from 
eravimetry, is included in this section. Comparative analyses of this sort 
do much to minimize the speculative nature of crustal interpretations 
produced from reduced gravity data alone. 

1. 160 kJ Seismic Profile Data 

Figure 19 depicts the tracklines over which USNS BARTLETT (T-AGOR 13) 

made runs employing a 160 kJ seismic sparker system. Continuous seismic 
profiles were obtained by this relatively low-resolution system. These 


were compared with the profiles obtained from the author's CBA isoline 
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chart (Fig. 15). The measurement units for the seismic profiles are 
one-way sound pulse travel time in milliseconds. This is more realistic 
than conversion of the crdinate values to depths since to do so would 
require the inferrence cf a density structure so as to calculate appro- 
priate sound velocities. This would only degrade the accuracy of the 
profile. The postulated basement depths are sketched in dashed lines 
since they have been deduced solely from seismic profiles. 

On trackiine A'A (Fig. 20) the discontinuity shown at Station 
80 corresponds to the Palo Colorado-San Gregorio fault zone, with a 
probable dip-slip component of about 430 m. The seismic data, in agreement 
with the fault map of Greene (1973), indicates this fault zone extends to 
within a few meters of the sea floor. The maximum CBA value of 14 mgal at 
Station 80 probably signifies the minimum depth to the dense granitic 
basement rock due to faulting. 

The two vertical lines at Station 81 indicate probable basement 
faulting as interpreted from the seismic profile. The radical change in 
the CBA profile agrees with the seismic data, in which a 3 mgal/km gradient 
is exhibited in the fault zone area. The general downward trend of the 
proposed granitic basement in the seismic profile is to some extent par- 
alleled in the decreasing CBA values to the west of Station 80. Also, the 
lower depth to basement to the east of Station 80 is reflected in the 
corresponding lower CBA values. 

Trackline AA (Fig. 21) shows excellent agreement between pro- 
posed granitic basement depth and the plotted CBA values. The two bulges 
in the otherwise smooth CBA trend coincide with the basement faulting 


interpreted from the seismic profile. No major faulting is apparent in 
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Figure 20. Comparison of CBA and Seismic Profiles for 
Trackline A'A. (Upper profile interpretation 
by R. S. Andrews, NPS, oral commun: > 1973) 
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profile interpretation by H. G. 
Greene, USGS, oral commun., 1973) 
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the seismic data but the decrease in basement depth from west to east 
agrees with the northward trending slope of the basement from south of 
Santa Cruz to the previously mentioned onshore granitic outcropping. 

Figure 22 corresponds to track B'B in Figure 19. The prominent 
correlative feature in these drawings is again the Palo Colorado-San 
Gregorio fault zone at Station 71. The probable area of dip-slip motion 
agrees well with the maximum change in CBA values. The upthrown side of 
the fault zone shows a gravity high which most likely corresponds to the 
granitic basement being at a shallower depth than in the surrounding area. 
This in itself is a manifestation of the local seismic activity. The 
small irregularity in the seismic profile near Station 37 is probably a 
part of the Monterey Bay fault zone. No correlative feature is seen in 
the CBA profile. This is undoubtedly shee to the depth of the basement in 
this area; that is, the change in gravity due to the faulting was filtered 
out by the Plio-Pleistocene crustal structure above it. it is possible to 
interpret the small deviation in basement profile at Station 34 as being 
in coincidence with the CBA high for that area. 

Trackline C'C (Fig. 23) displays excellent correlation between 
the seismic profile data and the CBA vertical cross-section. First, the 
general trend and slope of MLne basement follows that of the CBA profile 
very closely. The Palo Colorado-San Gregorio fault zone corresponds to 
the high CBA gradient just west of Station 69. Again, the increase in the 
CBA further west agrees with the uplifted side of the fault zone, bringing 
denser basement granite closer to the crustal surface. The maximum CBA 
value of 21 mgal, at Station 5, is just to the south of Santa Cruz where 
the shallow basement is known to be sloping upwards to the north. This 


is also accurately depicted on the seismic profile to the east of Station 6. 
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yo 2 KJ DeismicwkcolLile Data 

Figure 24 Show segments of three different tracklines obtained 
from a seismic survey using a 23 kJ sparker system (Greene et al., 1973). 
This system resulted in about 1500 m of penetration with approximately 5 m 
resolution. Since the tracklines cross the author's survey area, com- 
parison can be made between Greene's seismic profile interpretations and 
the author's corresponding CBA profiles. 

The fault extending up to the sea floor in the bottom drawing of 
Figure 25 is at the eastern edge of the Palo Colorado-San Gregorio fault 
zone. The upper figure shows good correlation with the gravimetry as 
there is an abrupt change in the overall CBA trend at the fault location. 
Higher values of CBA and abnormal gradients coincide with it. The small 
subsurface fault depicted on the seismic profile at the northeast end of 
line FF' exhibits no characteristic modifications in the CBA trend. This 
eould be due to the fact that it is of small extent (Fig. 24) and, as 
can be seen in the seismic profile reproduction, does not appear to alter 
the granitic basement (continuous solid line in Fig. 25). Thus, the 
lack of any dip-slip motion precludes the movement of any of the base- 
ment rock closer to the surface where a break in the CBA trend might 
be noticeable» Further, it is observed that the CBA profile parallels 
the seismically determined granitic basement. 

Trackline GG' (Fig. 26) does not traverse the Palo Colorado- 
San Gregorio fault zone but does cut across the northwestern end of the 
Monterey Bay fault zone. Due to the intermediate~resolution characteristics 
of the sparker used, the granitic basement was not decypherable. However, 
the steep gradients of the CBA values (5 mgal/km) correspond with the high 


density of faults running through the area. As seen from the seismic 
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Figure 24. Geophysical Tracklines for 23 kJ Seismic Profile Data. 
(Fault zones after Greene ct al., 1973.) 
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profile, the three separate faults appear to displace some Holocene deposits 
at the seafloor but probably do not reach the crustal surface itself. The 
overall trend of the strata shown in the seismic profile lends credence to 
the assumption that the general basement slope is on the order of the in- 
clination exhibited by the CBA cross-section. 

Trackline HH’ (Fig. 27), as shown, extends to within 1 km of 
land to the east of Santa Cruz (Fig. 24). In fact, the high CBA values 
and rather horizontal structural characteristics indicated in the seismic 
profile for the northeastern 5 km of the track correspond to the northward 
trending upslope of the granitic basement from south of Santa Cruz (Cronyn, 
1973). Again, the CBA gradient correlates well with the one fault zone 
depicted on the seismic profile. That no abrupt increase in CBA value 
is evident is probably due to the fact that the fault does not reach up 
to depths shallower than that of Pliocene age. Also, the general trend of 
the CBA profile parallels the predominate slope of the subsurface strata, 
which in turn probably reflects the general basement slope. 

3. Well Core Data 

Ross and Brabb (1973) have tied in known well depths to basement 
rock and derived a map showing depths to the granitic basement for most 
of the Monterey Bay area. OS chodeh highly speculative in nature, it never- 
theless gives a starting point in attempting to portray the true basement 
depths. Line EE’ (Fig. 28) connects two drill sites and extends west- 
ward into the oceanic portion of the author's survey area. Line DD' was 
picked so as to cross perpendicular to EE’ and intersect one of the test 
wells. The Humble Oil and Refining Company test well reached granitic 
basement rock at approximately 228 wn below MSL, while at the Monterey Bay 


Oil Company site the basement lies at slightly greater than 1620 nm. 
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Figure 27. Comparison of CBA and Seismic Profiles for 
Trackline HH'. (Lower profile after Greene 
etal 197 oe) 
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Figure 26. Geographical Location of Profiles for CBA 
Comparison with Well Core Data. 
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Figure 29 shows the excellent correlation which exists between the 
proposed basement location and the corresponding CBA trend. Based on this 
evidence, verification of Ross and Brabb's predictions of basement depth 
along line DD' seems legitimate. Although no fault zones are depicted in 
Ross and Brabb's interpretation of the basement profile, some undoubtedly 
do exist as DD' cuts across part of the Monterey Bay fault zone. Further- 
more, the northernmost end of DD' stops at the theoretical extension of 
the King City fault. The steep gradient (4 mgal/km) exhibited in the 
CBA values is usually indicative of subsurface deformation. 

The western end of line EE' (Fig. 28) lies at the eastern edge 
of the Palo Colorado-San Gregorio fault Zone; the other end is at the 
intersection of the forementioned King City fault and the conjectured 
Santa Cruz fault. Figure 30 illustrates the comparison between Ross and 
Brabb's proposed basement depths and the author's corresponding CBA cross- 
section profile for EE’. Good agreement exists for the area removed from 
the fault zones. The lack of agreement at the western end of line EE' 
is not significant due to its remoteness from the Humble Well, coupled 
with the fact that near the Palo Colorado-San Gregorio fault zone, in- 
congruous CBA profiles have been the rule. Due to the fact that faulting 
is in evidence near point E', the difference between mapped basement depth 
and the CBA profile there appears realistic. Previous figures have shown 
sharp increases in CBA values over fault zone areas and this seems to be 
the case near the Monterey Bay Oil weil. 

4. Sea Surface Gravity Data 
Sea surface gravity measurements, when corrected for horizontal 


and vertical ship accelerations, cross-coupling effects, instrument drift, 
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earth tides, earth curvature, and the Eotvos effect (due to changes in 
centripetal acceleration), and with theoretical gravity removed, give 
the free-air anomaly. 
The BARTLETT survey which produced the 160 kJ seismic data 
previously discussed, also collected sea surface gravimetric data. 
Gravity values for tracklines BB' and CC' (Fig. 19) are plotted for 
comparative purposes. As is seen in Figure 3l, the sea surface data con- 
tains large fluctuations in free-air anomaly values. These can be attributed 
to the high sea state which was encountered during this portion of the cruise. 
The Palo Colorado-San Gregorio fault zone shows up well on the 
BARTLETT plott and corresponds to the more accurate MFAA profile of the 
ACANIA data. Between the 5 and 10 km marks the sea surface data is totally 
unreliable; there is no correlation at all with the seafloor values. The 
Monterey Bay fault zone is clearly indicated on both profiles. Overall 
sea surface FAA values were two to three times greater than those for the 
bottom survey. Acceleration reinforcement between ship and meter oscilla- 
tions during the surface gravity measurements probably accounts for this. 
Figure 32 shows better correlation between sea surface and seafloor 
measurements but again some of the values obtained from the BARTLETT survey 
are two to three times greater than the corresponding ACANIA data. The 
two profiles do illustrate that the horizontal location of the maximum and 
minimum points are essentially coincidental. The horizontal differences 
observed at 7 and 15 km probably are due to unreliable navigation during 
the poor weather on the BARTLETT cruise. The large-scale geological 
features (Palo Colorado-San Gregorio and Monterey Bay fault zones) are 


explicitly portrayed on both profiles. The sea surface FAA is in better 
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Figure 32. Comparison of Sea Surface and Seafloor Gravity Data 
for irackiline Ce’. 
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agreement with the corresponding seafloor MFAA for trackline CC' than BB' 
probably because the former trackline cuts across the crustal deformations 
at a greater angle than the latter. 

5. Sea Surface Magnetic Data 

On the same cruise that the forementioned data was obtained, 
residual magnetism was also measured. The regional magnetic trend was 
ascertained from the International Geomagnetic Reference Field and the 
differences between that and the measured quantities were calculated. 
Since all information reduced to negative values, it can be concluded that 
the area surveyed by the BARTLETT either had an absence of magnetic materials 
in the crust or reversed polarity was in evidence (Maxwell, 1971). 

Tracks BB' and CC' (Fig. 19) were used in comparing the proposed 
basement configurations interpreted from the 160 kJ seismic data (referred 
to the corresponding CBA profiles) and the magnetic data. Approximately 
4 km east of point B (Fig. 33) excellent correlation exists between the 
crustal assymetry of the Palo Colorado-San Gregorio fault zone and the 
jump in regional magnetism. The gradient of approximately 50 y/km is 
indicative of subsurface crustal faulting. Eastward, the magnetic values 
are fairly constant and apparently do not reflect the small fault zone at 
11 km because of its depth and minimal vertical basement displacement. 

Figure 34 shows the same high order of correlation between the 
basement profile (verified to a certain extent by the CBA profile) and 
the residual magnetism cross-section plotted for trackline CC'. The 
higher values of magnetism near the westward end of CC’ probably reflect 
the approach to the edge of the continental shelf where negative magnetic 
values may not be so predominate (Jakosky, 1957). In this case there 


exists a magnetic gradient of approximately 70 y/km at the Palo Colorado- 
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San Gregorio fault zone. Again, this increase from the values recorded 
across track BB' can undoubtedly be attributed to the fact that line CC' 
crosses the Palo Colorado-San Gregorio fault zone almost perpendicularly. 
The absence of any basement faulting is represented as a fairly constant 


magnetic trace eastward from the large fault zone. 
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VI. CONCLUSIONS 


There are two major areas of faulting in the area investigated in 
this report. The most prominent is the Palo Colorado-San Gregorio fault 
zone: a narrow, northwest-southeast trending zone that joins the San 
Gregorio fault onland at Ano Nuevo Point. Through interpretation of the 
author's data, there is support for the theory that a fault of this 
magnitude can be traced from horizontal and vertical plots of CBA values. 
In fact, it appears that CBA irregularities will almost always become 
manifest if the granitic structure of the basement is vertically dis- 
placed (Fig. 17). 

The other zone, the Monterey Bay fault zone, parallels the Salinas 
Valley and the Sierra de Salinas. It is a wide belt of faults that 
crosses the southeastern portion of the author's survey area and closely 
approaches, but does not appear to cross, the Palo Colorado-San Gregorio - 
fault zone. CBA profile analysis aids in verification of seismic data 
in broad fault regions such as this but CBA data alone is of little use 
in this case except to trace likely areas of faulting. 

Fault displacement along the Palo Colorado-San Gregorio fault zone 
is similar to that on the San Andreas fault, that is vertical dip-slip 
and horizontal strike-slip motion, with rocks on the seaward side being 
displaced to the north. The prediction that an earthquake of magnitude 
7.2 to 7.9 could occur on the Palo Colorado-San Gregorio fault zone 
(Greene et al., 1973) lends much importance to corroborations of seismic 


profile data through results of gravimetric surveys. 


104 





VII. FUTURE WORK 


Because gravimetric analysis does not provide a unique model of the 
earth's crustal structure, corresponding interpretation of seismic and 
magnetic data is called for. It is recommended that seafloor gravimetry 
be extended south along the Palo Colorado~-San Gregorio fault zone to the 
maximum attainable depths. This should then be tied in with sea surface 
gravity and magnetic values to verify or modify existing seismic records. 

Since excessive depths and/or slopes in the vicinity of the Monterey 
Submarine Canyon preclude the possibility of bottom gravimetry in that 
area, sea surface gravity data and corresponding magnetic information 
should be correlated with previously interpreted seismic profiles. 

The ultimate goal of geophysical investigations in this marine 
environment is to produce a continuous picture of the off-shore CBA for 
the entire Monterey Bay area. Incorporation of all available data and 
maximum utilization of seafloor gravimetry should, in the near future, 


yield this objective. 
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